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Analytical expressions for the exchange interaction between J-multiplets of interacting metallic 
centers are derived on the basis of a complete electronic model which includes the intrasite relativistic 
effects. A common belief that this interaction can be approximated by an isotropic form oc Ji • J 2 
(or oc Ji • S 2 in the case of interaction with an isotropic spin) is found to be ungrounded. It is also 
shown that the often used “1 /U approximation” for the description of the kinetic contribution of 
the exchange interaction is not valid in the case of J-multiplets. The developed theory can be used 
for microscopic description of exchange interaction in materials containing lanthanides, actinides 
and some transition metal ions. 
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I. INTRODUCTION 

Strong magnetic anisotropy induced by spin-orbit cou¬ 
pling on the metal sites is a key ingredient for a num¬ 
ber of intriguing properties of magnetic materials, such 
as single-molecule magnet behavior mm, magnetic mul¬ 
tipole ordering 3], and various exotic electronic phases 
mm. If the spin-orbit coupling exceeds the crystal-field 
splitting of the ground term LS on the metal site, the 
latter acquires unquenched orbital momentum L and the 
low-lying spectrum is well described as crystal-field split 
eigenstates of the total angular momentum J = L + S, 
where S is the spin of the metallic term. This situa¬ 
tion takes place in lanthanides [6], actinides [3 and some 
transition metal ions in a cubic symmetry environment 

[2 IS]. 

The exchange interaction between such split-J crystal- 
field levels (or groups of levels) is significantly more com¬ 
plicated than the exchange interaction between pure spin 
terms (L = 0) described by the Heisenberg Hamiltonian 
J Si • S2. For the weak spin-orbit coupling, the dis¬ 
crepancy of the exchange Hamiltonian from the isotropic 
form was first pointed out by Stevens [9]. Later, the 
anisotropic exchange interaction was extensively devel¬ 
oped by Moriya m based on the Anderson’s microscopic 
approach mm- In the case of the strong spin-orbit cou¬ 
pling, the exchange interaction including the higher order 
terms of J was also phenomenologically treated since long 
time ago mm- The microscopic description was ad¬ 
dressed for the first time by Elliott and Thorpe [15] for 
uranium oxides, and by Hartmann-Boutron m for tran¬ 
sition metal compounds on the basis of simplified analysis 
based on so-called 1/U approximation (U is the average 
electron promotion energy between the sites). Recently, 
within the same approximation, the microscopic deriva¬ 
tion of the exchange Hamiltonian between J-multiplets 
was completed by Santini et al [3]. 

Despite this early evidence of complexity of exchange 
interaction between metal ions with unquenched orbital 
moments, it was repeatedly conjectured that the ex¬ 
change interaction between fully degenerate J-shells, in¬ 
volving (2Ji + l) and (2J2 + I) angular momentum eigen¬ 


states on the first and the second magnetic centers, re¬ 
spectively, is described by an isotropic exchange Hamil¬ 
tonian written in terms of J i momenta: 

^Heis = JJi • J 2 • (1) 

Contrary to Heisenberg Hamiltonian for isotropic spins, 
there is no a priori justification for the Hamiltonian 0 - 
Nevertheless, this form is often used for the description of 
interaction between lanthanides or actinides (or a similar 
form, oc Ji • S2, in the case of their interaction with an 
isotropic spin), especially, in the last years [T7H3Q] . One 
of the reasons that the simple bilinear form has been often 
used is that the large numbers of the phenomenological 
exchange parameters cannot be easily determined. 

It is not clear, however, how important are “non- 
Heisenberg” terms in the actual J-J coupling, nor is the 
1/U approximation a priori justified for metal ions with 
unquenched orbital momenta. Both these questions can 
only be answered after a more complete derivation of ex¬ 
change interaction between J multiplets on the basis of 
a reliable microscopic model. Besides, a microscopic de¬ 
scription of J-J ( J-S ) exchange interaction is desirable 
due to a very large number of phenomenological parame¬ 
ters, in contrast to weakly-anisotropic spin systems con¬ 
taining only a few of them [ED 02]. Given that many 
microscopic electronic parameters describing individual 
magnetic centers and their interaction can be accurately 
derived via density functional theory [33] or ab initio cal¬ 
culations m, a microscopically derived Hamiltonian for 
J multiplets can become a powerful tool for the inves¬ 
tigation of exchange interaction in materials containing 
lanthanides, actinides and transition metal ions with un¬ 
quenched orbital momentum. To this end the electronic 
Hamiltonians only need to be downfolded on the reduced 
manifold of low-lying states at the corresponding metal 
ions. 

In this work, we derive analytically the exchange 
Hamiltonians for interacting J-multiplets and for inter¬ 
acting J-multiplet and isotropic spin, starting from a mi¬ 
croscopic electronic Hamiltonian including the relativis¬ 
tic interactions on the metal sites. The obtained ex¬ 
change parameters are expressed via electronic matrix 
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elements which can be derived from electronic structure 
calculations. The structure of the exchange Hamilto¬ 
nian is discussed and the result is applied for the anal¬ 
ysis of some systems with different geometries. Com- 

g rison with the predictions given by the Hamiltonian 
and the simplified treatment on the basis of 1/U ap¬ 
proximation shows that both of them are not suitable 
approaches to describe the exchange interaction of ions 
with unquenched orbital momentum. Finally, the rel¬ 
ative contributions to the kinetic exchange interaction 
from intermediate states is analyzed. 


II. MICROSCOPIC DESCRIPTION OF 
INTERSITE INTERACTION 

We derive the expression for the interaction between 
metal ions with unquenched orbital moments. The 
derivation is based on a complete electronic Hamilto¬ 
nian, including all intrasite relativistic effects, and em¬ 
ploys adequate approximations. The multipolar intersite 
interactions of electromagnetic type, such as the elec¬ 
tric quadrupolar and magnetic dipolar interactions, have 
been described elsewhere [3, 14, 35, 36 and are not con¬ 
sidered here. Their effect can be taken into account as 
additive contribution to the exchange parameters. 


A. Electronic multiplets on sites 


The nonrelativistic electronic state of an ion with par¬ 
tially filled nl shell (n is the main quantum number and 
l is the one-electron orbital angular momentum) corre¬ 
sponds to an LS'-term characterized by the total orbital 
L and spin S angular momenta [37|- The eigenfunc¬ 
tions {| (XlsLMlSMs}} are described by orbital and spin 
quantum numbers, L and S', and the projections of L and 
S on a given axis z, Ml and Ms, respectively; qlls indi¬ 
cates the other quantum numbers. The (2 L + 1)(2S +1)- 
fold degenerate term is further split by the spin-orbit 
interaction into J-multiplets which are eigenstates of the 
total angular momentum J = L + S. The corresponding 
eigenfunctions {| ajJM)} are characterized by quantum 
numbers of the total angular momentum J and its projec¬ 
tion M = — J, — J+1, • • • J (aj stands for other quantum 
numbers). 

In general, the spin-orbit interaction mixes multiplets 
with the same J belonging to different LS terms (the 
so-called J-J mixing [36]). In the case when this mixing 
can be neglected, each J-multiplet is attributed to one 
LS -term and the corresponding wave functions become 
of the form (aj = {qlls, L, S)): 


\otjJM) = Y, I cxlsLM l SMs)C j l % lSMl , (2) 
M L M s 


where Clm l sm l * s the Clebsch-Gordan coefficient (A2) 
[38] . This is a good approximation, in particular, tor 


the ground J-multiplet of trivalent ions from the late 
lanthanides series, Ln 3+ . 

When the metal ions are embedded in complexes or 
crystals, their electronic structure is modified due to co¬ 
valent and electrostatic interaction of the magnetic nl or¬ 
bitals with the environment. In the case of /-metals the 
magnetic orbitals are usually strongly localized and the 
effect of the surrounding is relatively weak. For exam¬ 
ple, in the case of lanthanide, the intraionic bielectronic 
interaction leading to atomic terms separation is ca 5-7 
eV and the spin-orbit splitting is ca 1 eV for lanthanide 
ions, thus exceeding several times the the crystal-field 
splitting, which is usually of the order of 0.1 eV m- 
In this situation the low-energy electronic states are well 
approximated as crystal-field split atomic J-multiplets. 
Due to the weak hybridization of the 4/ orbitals with the 
surrounding, the Wannier functions of the corresponding 
magnetic orbitals practically coincide with the atomic 4/ 
orbitals. Similar holds true for actinide ions although 5/ 
orbitals are more delocalized than 4/ orbitals. 

On the other hand, the effect of the hybridization of d 
orbitals with the ligand orbitals is usually much stronger 
than in lanthanides and actinides resulting in a crystal- 
field splitting which often overcomes the atomic LS-term 
splitting. Therefore, the orbital angular momentum is 
generally not a good quantum number for nondegenerate 
ground state of embedded transition metal ions. More¬ 
over this orbital angular momentum is quenched as a rule, 
(L) = 0, in most compounds. The exception is the cubic 
environment, in which the d orbitals split into doubly de¬ 
generate e and triply degenerate £2 levels. When the £2 
orbital are partially filled, the electronic state is charac¬ 
terized by the nonzero fictitious orbital angular momen¬ 
tum L = 1, which couples to the total spin of the site via 
spin-orbit coupling and gives molecular multiplets char¬ 
acterized by fictitious total angular momentum J = L + S 
(see Ref. [8]for details). 


B. Intersite interaction 


The electronic Hamiltonian H for electrons localized at 
two sites can be divided into the intrasite contributions 
J£q = 1,2), the intersite bielectronic iJjb and electron 
transfer H t parts: 

H= Y H'o + K+Ht. (3) 

i= 1,2 


The intrasite Hamiltonian for site i, iJg, contains all ef¬ 
fects discussed in Sec. Ella] such as the non-relativistic 
atomic terms, the spin-orbit term and other relativistic 
corrections, and the crystal-field. The eigenstate of Hq 
is determined by the number of electrons Ni in magnetic 
orbitals and crystal-field level p, \iN i: p). i/b consists of 
intersite Coulomb interaction iJcoui and direct exchange 
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(multipole) part i^DE- 

#bi = ^Coul + ^DE, 

Pcoul — ^ ^ U fllm<jft J 2ncr' •> 


( 4 ) 

( 5 ) 


^de = — 


Vmm'n'n^lma^ria' ^ 2 m'a' ^2n 'a •> (6) 


mnm' n' era' 


where ra, n indicate the projection of the orbital angu¬ 
lar momentum l i: a is the projection of the electron spin 
momentum, cj ma (Q mcr ) is the electron creation (annihi¬ 
lation) operator in spin-orbital (ra, cr) of site i (= 1,2), 
— ^ , U' is the intersite electron repulsion, 


n. 


= c. 


'ima^ irncr ’ 


Vmm'n'n 


is the inter site exchange integral, 

= J dr 1 dr 2 ^ m (r 1 )ip^ m ,(r 2 ) 


Vr. 


mm'n'n 


x vflri - r 2 |)'02n(ri)V’in'(r2), 


( 7 ) 


'ijjim is the Wannier function at site i and component m, 
and v(|ri — r 2 1) is the two body interaction. Note that 
H[ a does not change the number of the electrons on each 
sites. The transfer Hamiltonian is written as 

Ht = ^mm'^lma^jm'ai ( 8 ) 

i^j mm'a 


where is the electron transfer parameter between 

orbitals m of site i and m' of j. 

As discussed by Anderson [lTJ[T2], the direct exchange 
parameter 0 and the transfer parameters are finite due 
to the delocalization of the magnetic orbital on the 
ligand between metal sites (or the other magnetic site). 
The electron transfer parameter between metal sites is 
at least several times smaller than the intrasite electron 
repulsion. Thus, the electronic Hamiltonian H ([ 3 ]) can 
be divided into zeroth order Hamiltonian 


where r indicates the eigenstate of the system and Cp q ]^ 2 
is a coefficient. In the derivation of the effective Hamil¬ 
tonian, we consider the truncated vector space B of the 
electron configurations which include up to one electron 
transfer with respect to the numbers of the electrons in 
the ground electron configurations. Hereafter, Ni is used 
as the number of electrons on site i in the ground electron 
configurations. For simplicity, the numbers of electrons 
in the ground configurations will not be written explic¬ 
itly, and the configurations 8 with Ni — 1 and IV 2 + 1 
(Ni H- 1 and N 2 — 1) electrons on sites 1 and 2 are ex¬ 
pressed by the type of the virtual electron transfer, 1 —2 
(2 —)> 1). Therefore, B is defined as follows: 

B = {| r )> |1 2, s), 12 1,8')}. (13) 


The eigenenergies of the states |r), |i j,s) are denoted 

as E® and E l s ~^i , respectively. 

The exponent of the unitary operator e s is given as 



P^P t P r Q 
P° - E^ j 


p y °p t prA 

EO-E*r j ) 


(14) 


where P° and Pj^- 7 are the projection operators, 

P? = \r)(r\, (15) 

PV 3 = M h s)(i ->• j, s\. (16) 


The exponent S is chosen to fulfill the condition 

[S,H 0 + HDE]=H t (17) 

within the space B. The effective Hamiltonian ( fIo| ) 
within Bo = {\r)} is obtained as 


H eS = H 0 + H DE -~[S,H t ], 


(18) 


2 

ff 0 = 5Z#0 + tfcoui (9) 

2 = 1 

and small terms Pde and H t . For localized magnetic 
electrons, the latter can be treated in the first and the 
second order of perturbation theory, respectively mm- 
This is done here via a unitary transformation, 

H eff = e~ § He § , (10) 

removing H t from the initial Hamiltonian. Neglecting 

the terms higher than second order after the transfer pa¬ 
rameters, we obtain the effective Hamiltonian acting on 
the ground J-multiplets on sites, 

P eff = P 0 + P. (11) 

In the unperturbed Hamiltonian ([9|, there are no terms 
which vary the numbers of the electrons. Then the elec¬ 
tronic states can be written as follows: 

\N u N 2 ,r) = '52\lN 1 ,p)\2N 2 ,q)C&?*, (12) 

p,q 


up to second order after H t . The second and the third 
terms in Eq. (18) correspond to H defined above, 


H = P D e + Pke, 


*- = *EEE 

r s 


E° - Ei^ 


+ H.c. 


(19) 

( 20 ) 


Note that the terms and H t P^H t P’ s ^ 7 do 

not enter here because they map the states outside the 
domain Bo. 

Neglecting the effects of the crystal-field splitting in 
the denominator of Pke> which is a reasonable approxi¬ 
mation for our systems, the eigenstates r, s reduce to the 
sets of the J-multiplet states: 


r {ol\J\M\,ol 2 J 2 M 2 ), (21) 

( 22 ) 

where, Ji is the total angular momentum with the ground 
electron configuration, and J, J r are the total angular mo¬ 
menta for intermediate states arising from the transfer 
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of one electron between the sites. The kinetic exchange 
Hamiltonian becomes 


a^ = \E E E 

aiJi,ajJj ajJ,a'jJ 1 

+ H.c., 




E P° 


t r a j n t r oiiJi,aj Jj 

E° — 

^OiiJi,ajJj - Cj ajJ,a'jJ' 


where the projection operators are 


pO _ p 

1 Oi-i Ji ,CXj J j 1 " 


,Ni 




r cxj 


and P n N 


iociJi jotjJji 

yNi-lpNj+l 
iajJ jct'jJ n 

iajJ is the projection operator on site i, 


(23) 

(24) 

(25) 


jj= E \iNajJM){iNajJM\. (26) 

M— — J 

In the space of the ground J-multiplets on sites, Bj = 
J 2 M 2 ) : — Ji < Mi < J^}, the kinetic exchange 
Hamiltonian reduces to 


HtP^i / j,H t 

1 OLjJ,OLjJ' 1 

pO _ 

OLjJ.Oi'J' Ji,Jj oljJ,ol'jJ' 


^ke = E E 


(27) 


Here, oti of the ground J-multiplet is not written for the 
sake of simplicity and P^.j. a .j. is omitted because it is 
the unit operator within Bj. 

Substituting Eq. © into Eq. ([27]), we obtain 

HkE = E E E E E “Cm'Cn 

i^j ocjJ ol'jJ' mna m'n'< j' 

(c t ,)fc- / P N J + 1 c f 

\jim<j r %olj J Cmcr J \pjm cr r ja'jJ'^jn'cr' J 


Ua + A E%~} + A E Nj , + j, 

L J lOCjJ 1 JOijJ' 


,(28) 


where Uij is the smallest promotion energy for the elec¬ 
tron transfer from site i to site j, and A E^ ajJ is the ex¬ 
citation energy from the ground intermediate state with 
N electrons. In the derivation of Eq. (28), we have 
not used the approximate form ([2]) for the J-multiplets 
in both the ground and the virtual states. The quan¬ 
tum number J for the virtual states fulfills the condi¬ 
tion \Ji — li — l/2\ < J < Ji + li + 1/2. Although the 
crystal-field splitting is neglected, the multiplet struc¬ 
tures of sites are completely retained in Eq. (28). Hence, 


important effects such as the Goodenough’s mechanism 
40 are included in the kinetic exchange Hamiltonian. 


III. 


EXCHANGE HAMILTONIAN IN 
J-REPRESENTATION 


The Hamiltonian H (19) is transformed into tensor 


form with the use of irreducible (double) tensor technique 
[41]. method of equivalent operator |8], and the form of 
J-multiplet state ©: 


H ^ ) iJkqk'i 

kqk' q' 


qUEoUhl 


(29) 


Here, O k ( Ji) and J 2 ) are Stevens operators [42] (see 

also Appendix |A 3| ) whose ranks k and k' have to obey the 
relation k + k' = even due to the invariance of the Hamil¬ 
tonian with respect to time inversion [8] (see Appendix 
[Cl ), q and q' are component, and O^(Ji) and O k ,{J 2 ) are 
the scalars obtained by replacing J? and J^ in O k (Ji) 
with eigenvalues JpJiE 1) and J^, respectively {i = 1,2). 
The exchange coupling constant Jkqk'q' is a sum of the 
direct exchange J B q E , q , and the kinetic J E q E , qf contribu¬ 
tions [431: 


Jkqk'q' ~ Jkqk'q' + Jkqk'q' * (30) 

The advantages of using the exchange Hamiltonian in 
the tensorial form is that with Eq. (29) it is easier (i) to 


obtain physical insight on the exchange interaction and 
(ii) to combine it with other terms such as crystal-field 
and Zeeman interaction included in Hq. The latter can 
be treated at ab initio level [321 [34] . 


A. Direct exchange interaction 

The outline of the derivation for the direct exchange 
part in H is given here, whereas details of the calculations 
are given in Appendix m The direct product of the 
double tensors cj mcr Q mcr / appearing in iJoE ^ is reduced 
as follows: 




(- 1 ) 


b+^+ b +cr 7 


aab /3 




v C iaa 
X °Lml 


pib/3 


(31) 


where c in(r 
tensor (see Appendix 


= (—l)h +™+ 2 +(j/ Cj-n-cr' (A16) is a double 
the curly bracket {cj (8) Ci}l 


A 2), the curly bracket {cj (g) 
indicates the irreducible operator of ranks (a, b ) and com¬ 
ponents (a, /3) constructed from the product of two ten¬ 
sors, where the superscripts and subscripts are the orbital 
and spin parts, respectively. The irreducible tensor op¬ 
erator is replaced by the total angular momentum opera¬ 
tor using the method of equivalent operator for double 
tensor ( |A26l ): 


kq 


T)' 

aab /3 u abk 


quE 

o° k (JiY 


(32) 


Here, T> l abk is a tensor with three indices a, 6, k ( |B2|). N ote 
that when the method of equivalent operator ( |A26 ) is 
used, the form of | J M) pi ) is assumed. 

Substituting Eqs. (plFf (32) into Hde (p| 5 we obtain 


the tensor form of the direct exchange Hamiltonian. The 
exchange parameter J BE is obtained as 


Ji 


DE 
kqk' q' 


' Vkqk'q r 'E > abk'D‘i'bk f i (33) 

aa'b 
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where 


yjaa'b _ \ A \ A \ ^ (_-\\li+l 2 +n+n'— fiy 
V kqk'q' ~ /.A 2-^ 2-^ ^ v mm'n'r 

mn m'n' aa'(3 


X c: 


C kq s~ia a* 


rik'q' 

1 n r rv 1 


Iimli—n^aab(3^l2m'l2—n'^a'a'b—f3' 


B. Kinetic exchange interaction 


(34) 


The derivation of the kinetic exchange parameter J KE 
is similar to that of J7 DE . As in the previous case, only 
the outline of the derivation is given here and the de¬ 
tails can be found in Appendix |B 2| In comparison with 
the direct exchange, the derivation of J KE is cumber¬ 
some because of the projection operator P^jj appearing 
in Eq. ( [28] ) . The latter is reducible within the product 
group S' 0(3) (8 ) SU (2) where the creation $ and annihi¬ 
lation c ( |A16 ) operators are irreducible. Therefore, we 


first transform PP T into the sum of the irreducible dou- 

LUCJ J 


as ct 


ble tensors T > A jaaa , a , (Bll), and then the product such 


ima iajaaa'a' ' incr ‘ 


,Ci na > is reduced. The irreducible tensor 


operator {cj (8) {Pia'jaa ® c i}d}el * s replaced by the total 
angular momentum using the method of the equiv¬ 
alent operator (A26), and finally we obtain the kinetic 


exchange parameter: 


‘TfcgfeV - E E E 


X ^kqk'q'^ajJfxkSh'jJ'fx'k’ 


fxx' OLjJ Oi'jJ' 

EEE 

fxx' ajJ ol'J' 


u 12 




\fxx r 


\ & 
' n' ^ cm 


A E: 


F 2 


'N% -h 1 


\t X £} kqk' q/'d aj jf x k J a'jJ'fx'k' 

U 2 i + A E[ 


iVi+1 

OLjJ 




(35) 


where 


{t X nlqk'q' = {-l) h ~ h - f+q ' E E E 


+12 .21 

o 1 v m 


I 1 —I 2 — f ~\~q 

mn m'n' 

r-ixf y-ixf, r^x'f,' 

A ^hnkq^fcfhm^h-n'k'q'^f-(f>l 2 -m'' V°°/ 


The range of variation of indices of the tensors described 
above, as we ll as in the subscripts of T and Q is specified 
in Appendix B 2 As in the case of the derivation of J EE , 
Eq. ([ 2 ]) is used. Because of this assumption, the quantum 
numbers of the intermediate multiplets contributing to 
J KE obey the relations: \Li — li\ < L < Li + l i: |S^ — 
1/21 < S < Si + 1/2, and |Ji——1/2| < J < Ji+h + 1/2, 
where L, S', J without subscript i refer to the intermediate 
states l44l. 


the rank k and k' in Eq. (29) is only determined by the 
electronic state of site 1 and 2, respectively [45] : 


tax = min[2/i + 1,2J/. 


(37) 


Thus the maximum rank for /-electron system, f N , is 7 
for N = 2-4 and 7-13, & max = 5 for N = 1,5 and /c max 
0 for N = 6 [46] . 

On the other hand, the range of q is determined by the 
nonzero parameters describing the intersite interactions, 
Vmm'n'n and If A max is the maximal difference 

of the indices corresponding to one site in the above pa¬ 
rameters (m — n for site 1 and m' — n' for site 2) then 
the upper bound for q ( q ') is 


g m ax = min [A 
Note that terms with 


: T 1, ^max] • 


(38) 


^max 


will also be present in the 
Hamiltonian (29) due to the time reversal symmetry, im¬ 
plying the following range: -g max < q < g max . 

The effective Hamiltonian (29) is further divided into 
the exchange H ex and the zero-field splitting parts. The 
latter is defined as comprising terms with either k = 0 or 
k' = 0. 


D. Decomposition of H 


Knowledge of the domain of rank k (37) in the general 


exchange Hamiltonian (29) allows us to calculate Jkqk'q' 
by using the orthogonal] 


ity of the Stevens operators: 


yJkqk'q' — ( 1) 


q+q 


n 


n Jij 2 
X Tr — 


kk' /O 1 / 2 J 2 

Ji/c0 o J 2 J 2 /c / 0 


(39) 


where = y/2j + 1 and Ujjt = YljUjf, the trace (Tr) is 
taken over the ground J-multiplets, and 


Qkqk,q '- otwoUkY 


(40) 


The form (39) for the exchange parameters offer some 
advantages for practical calculations. H enters Eq. 
(39) in the form of the numerical matrices in the ba¬ 


sis of the products of multiplet wave functions on sites, 
|JiMi, J 2 M 2 ). The exchange parameters obtained by 
Eqs. (33), (35) and those calculated by the projection 


(39) have been compared with each other for some test 


examples. 


IV. EXCHANGE INTERACTION BETWEEN 
J-MULTIPLET AND ISOTROPIC SPIN 


C. Structure of the exchange Hamiltonian 

The domains of variation of k and q characterize the 
structure of the Hamiltonian (29). The upper bound for 


When the orbital angular momentum is zero in the 
ground state term of one of the sites, the low-energy 
states of this site are characterized by the correspond¬ 
ing spin S. This situation is encountered in mixed 
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FIG. 1. (Color online) (a) Linearly bridged Dy dimer with 
one ligand atom (L). (b) Kinetic exchange interaction be¬ 
tween 4/( 5z 2_ 3r 2) z orbitals, (c) Calculated exchange spec¬ 
trum with full H ex (43) and with its different contributions 
(U = U 12 = U 21 = 5 eV). The excitation energies of the in¬ 
termediate state s entering Eq. (35) have been calculated ab 
initio (Sec. VA). 


lanthanide-transition metal and lanthanide-radical com¬ 
plexes men e]. The exchange Hamiltonian between a 
J-multiplet and an isotropic spin is obtained in a similar 
way as Eq. (29): 


H — ^ ^ \JkqQ0 

kq 


Q£(JiA 

o° k (J 1 ) 


$kqlq 


kqq' 


O q k(Jl)S2q' 

' OIWS2 • 


(41) 


V. EXAMPLES 


We further consider some typical examples of J-J and 
J-P exchange interaction. Since the kinetic exchange in¬ 
teraction is usually much stronger than the direct ex¬ 
change interaction puna, we only take into account the 
former. In order to include the multiplet structure of 
intermediate states in Eq. ( |35| ), first we calculate ab ini¬ 
tio the excitation energies of the virtual electron-transfer 
states. The calculated exchange levels are compared with 
those arising from the bilinear form 0 and correspond¬ 
ing to the 1/U approximation. 


A. Excitation energies of Dy + and Dy + ions 


The excitation energies A E aj j appearing in the de¬ 
nominator of the kinetic exchange Hamiltonian ( [28] ) are 
calculated ab initio. Since the effect of the crystal-field 
split ting in the intermediate states is negligible (Sec. 
IIA), we used the energy levels of the free Ln 2+ and Ln 4+ 


ions. In this work, we only calculated the energies for Ln 
= Dy that we will use in the following sections. Apart 
from the crystal-field splitting, there is totally symmetric 
electrostatic potential which only depends on the num¬ 
ber of electrons and shifts uniformly the J-multiplet en¬ 
ergies. This effect is absorbed in the minimum promo¬ 
tion energy Uij in Eq. (35). The energies are calcu¬ 


lated using the complete active space self-consistent field 
(CASSCF) and the restricted active space SCF state in¬ 
teraction (RASSI) methods with ANO-RCC QZP basis 
set [49]. With the CASSCF method, the PP-term en¬ 
ergies are obtained, while with the RASSI method, the 
spin-orbit (J-multiplet) energy levels are calculated. For 
the CASSCF calculations, all 4/ orbitals are included 
into the active space. The terms included in the RASSI 
mixing are 6 P, 6 P, 6 H for Dy 2+ and 7 P, 5 P, 5 P, three 
5 P>, two 5 P, three 5 G, two 5 P, two 5 /, 5 A, 5 P for Dy 4 +. 
The excitation energies are tabulated in Table [I] 


The expressions for exchange coupling constants are sim¬ 
ilar to Eqs. (33), (35) and are listed in Appendix |D| Be¬ 
cause of the lack of orbital degrees of freedom on site 2 
(I 2 = 0), the rank k' of the spin operator does not ex¬ 
ceed 1. Due to the time-reversal symmetry, k is even and 
odd for the first and the second terms in Eq. ( [41] ) , re¬ 
spectively. As in the previous case, the former (k r = 0) 
is zero-field splitting term and the latter (k r = 1) is the 
exchange interaction. 


There are several PP-terms which appear more than 
once, i.e., 5 P>, 5 P, 5 G, 5 P, 5 1 terms of Dy 4+ ion. These 
terms obtained by the CASSCF calculations are assigned 
to the symmetrized PP-states within the shell model, 

I l N ollsLSJM) , comparing the patterns of the ab ini¬ 
tio and the model spin-orbit splittings of each PP-term. 
The symmetrized PP-states are constructed by using the 
coefficient of fractional parentage [48]. In the basis of 
symmetrized states, the matrix element of the spin-orbit 
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TABLE I. Excitation energies of Dy 2+ and Dy 4+ (meV). The number in the parenthesis indicates quantum number aLS in the 
main text, where the enumeration follows Ref. 1481 


LS'-term 

Dy 2 + 

J 

A E& jj 

LS'-term 

Dy 4+ 

J 

A Ecx jj 

LS'-term 

Dy 4+ 

J 

A E a jj 

5 i 

8 

0.000 

7 F 

6 

0.000 

5 H{2) 

4 

8738.547 


7 

458.212 


5 

309.504 

5 I( 1) 

8 

10579.798 


6 

859.148 


4 

508.583 


7 

10778.918 


5 

1202.808 

5 D( 1) 

4 

13934.562 


6 

10917.180 


4 

1489.190 

5 D{ 2) 

4 

7624.320 


5 

11003.924 

5 D 

4 

5667.150 

5 D{ 3) 

4 

3298.723 


4 

11049.978 

5 f 

5 

2369.357 

5 F(1) 

5 

5838.268 

5 /(2) 

8 

5827.337 


4 

2655.740 


4 

5936.668 


7 

6035.414 

5 g 

6 

3412.346 

5 F( 2) 

5 

11206.930 


6 

6173.415 


5 

3756.005 


4 

11118.121 


5 

6242.469 


4 

4042.388 

5 G(1) 

6 

6925.862 


4 

6261.863 





5 

7221.721 

5 K 

9 

6516.513 





4 

7423.960 


8 

6766.583 




5 G(2) 

6 

4540.837 


7 

6935.080 





5 

4539.168 


6 

7051.111 





4 

4613.355 


5 

7130.959 




5 G(3) 

6 

13458.495 

5 l 

10 

4277.003 





5 

13323.410 


9 

4402.240 





4 

13191.767 


8 

4537.021 




5 H( 2) 

7 

8096.483 


7 

4674.613 





6 

8392.826 


6 

4809.212 





5 

8603.858 





Hamiltonian H so = £ 'U ' (C >0) is given by 
(l 4l+2 ~ N a L sLSJMj\H so \l 4l+2 ~ N a / LS L / S'J'Mj) 


= —5ll’5ss'5jj'5mjM'X 


n iJV /Z(ITT) 

y/2 



> (l N a LS LS\\{X <g> c}\\\l N a' LS LS) 


(42) 

for TV < 21 + 1. Here, the curly bracket with 3x3 ele¬ 
ments are the 9 j symbol ( |A5[ ). Therefore, the spin-orbit 
splitting is proportional to the reduced matrix element 
of operator { c i (8) c}{™ m 50. • 


B. Kinetic exchange through monoatomic bridge 

As a simple example, consider an exchange-coupled 
Dy 3+ dimer with axial bridging geometry (Fig. Ilk). 
The largest transfer parameter (t) is expected between 
hz 3 -3r 2 z i m = 0) orbitals because of their sigma bonding 
to the p z orbital of the bridging ligand atom (Fig. ER 
Then, according to the rule (38), q max = 1, while Eq. 


Hamiltonian iL ex , after expanding the Stevens operators 
in Eq. (29), is 


H ex = K w + Ji • 3 2 K (2) + K (2) Ji • J 2 , (43) 


k {1) = E 

k,k '=0 

a-«« = £ 

k,k' = 1 


(44) 

(45) 


(36) gives q = —q'. The resulting form of the exchange 


where k + k'=eve n. We can see that, even in this sim¬ 
plest case, Hex does not reduce to the isotropic form 
(1) because Ising (oc K^) and mixed Ising-Heisenberg 
(oc K ( 2 )) terms, both involving high powers of momen¬ 
tum projection operators of two sites. The parameters 
K ™,, are tabulated in Table [ h| When the eigenvalue 
of Ji z is large, the higher order terms are significantly en¬ 
hanced and contribute to the exchange interaction rather 
than the bilinear term. As a result the exchange spec¬ 
tra calculated with the full Hamiltonian (43) and with its 
Heisenberg-type part 0 show large discrepancy between 
them (Fig. [ljc)) . The discrepancy is also seen in their 
eigenstates. The difference between the exchange states 
of (43) with those of Hpeis 0 is compared by expanding 
the former by the latter. The solution of HHeis fc> r two 
















FIG. 2. (Color online) (a) Dy dimer bridged via anion (D 2 h core symmetry). Large (purple) and small (red) balls are Dy 
and N, respectively, (b), (c) Kinetic exchange interaction between 4/( 5z 2_ r 2) cc orbitals and 4/ x 3_ 3xy 2 orbitals, respectively, via 


(d) Kinetic exchange interaction between the 4 f xyz orbital of Dy and the LUMO of Nj . (e) Calculated 


the HOMO of N^~ 

exchange spectrum with full H ex and its first-rank contribution (U = U\ 2 — U 2 \ = 5 eV), and in the 1/U approximation for 
kinetic exchange pattern (b). (f) Calculated exchange spectrum with full H ex and its first-rank contribution ( U 2 ± — 3 eV), and 
in the 1/U approximation for kinetic exchange pattern (d). The calculations in (e) and (f) involved exchange parameters (35) 
with excitation energies of the intermediate states on Dy evaluated ab initio. U in the calculation within 1/U approximation 
was chosen to reproduce the width of the spectrum for full H ex . 


site system is given 

|Ji 2 M 12 )= ]T (46) 

Mi, M2 

where J 12 and M i2 are the total angular momentum for 
the dimer and its projection. The low-energy exchange 
states of H ex (43) are written in the basis of {| J 12 M 12 )} 
as follows: 


—0.278523, -E' 3 , 4 , 5,6 = —0.243015, respectively, in the 
units of t 2 /C/ 12 . The low-energy exchange states are not 
necessarily mainly contributed by the ground state of the 
antiferromagnetic Hueis, \J 12 ^ 12 ) = |0,0). Therefore, 
we conclude that the Heisenberg form of the interaction 
is not adequate to describe the exchange interaction be¬ 
tween J-multiplets. 


\^i,A lg ) « 0.574|0,0) + 0.77312, 0) + 0.190|4,0) 

- 0.160|6,0), 

\^2,A 1u ) « 0.846|1,0) + 0.50513,0) — 0.161|7,0), 

|T 3 ,4, E 1u ) w 0.35211, ±1) + 0.768|3, ±1) + 0.522|5, ±1) 
+ 0.107|7, ±1), 

|*5,6, E lg ) w 0.63612, ±1) + 0.71214, ±1) + 0.294|6, ±1). 

Here, the irreducible representation T of is used, 

and the states |\Eq,T) belong to the eigenvalues Ei^ = 


C. Kinetic exchange through biatomic bridge 

Consider the exchange interaction in the Dv 3+ dimer 
bridged by the (n = 2,3) anion (Fig. |2ja)) [I5T[ . 
In the case of n = 2, 4/ electrons of Dy 3+ ions would 
transfer between the metal sites via the highest occupied 
molecular orbital (HOMO) of (Figs. 2^b),(c)). The 
HOMO overlaps with the /( 5 Z ^- r ^) x (|m| = 1) and the 
fx 3 -3xy 2 ( \m\ =3) metal orbitals, the former interaction 
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A £/eV 


A £/eV 


(c) 



(d) 


(N 

1 ? 


0.6 

0.5 

0.4 

0.3 

0.2 

0.1 


■ . L .t. . . 1 .1. - 
2 3 4 

A E/eV 


5 6 


FIG. 3. The contributions from the intermediate states w corresponding to the excitation energy A E (= A E^jj + A E^, Jt ) 

for (a) the full exchange interaction H ex and (b) the exchange within the 1 /U approximation H'_{ U for Dy 3+ dimer bridged by 
N| _ and for (c) the full exchange interaction and (d) the exchange within the l/U approximation for the Dy-radical system 
(n = 3). The units of w are t 2 /U 21 for (a), (c) and t 2 /U for (b), (d). 


being dominant. Hence, we only consider the electron 
transfer between the orbitals with \m\ = 1 (Fig. Em 
For them A max = 2 and we obtain according to Eq. ( [38] ) 
(/max = 3. Then iL ex will include powers of Ji± (= Ji x ± 
iJiy) for each center up to third order. 

Figure [2je) shows the calculated exchange spectrum 
for full -f/gx, and its first-rank contribution, and for one 
single promotion energy U (1 /U approximation) [3l ITS]. 
Although the first-rank contribution is bilinear in J^ 7 , it 
is not isotropic and the corresponding spectrum does not 
resemble the pattern of levels of Heisenberg-type Hamil¬ 
tonian _ 0 - Also the spectrum is quite different when 
the l/U approximation is applied. This approximation 
neglects the splitting of the LS- terms which exceeds sev¬ 
eral times the minimal electron promotion energy. As a 
result, the relative contributions to the exchange inter¬ 
action from various intermediate states are significantly 
modified. In order to see the variation of the contribu¬ 
tions from the intermediate states to the kinetic exchange 
interaction, we divide the kinetic exchange Hamiltonian 


as follows: 

^ = EE h ex (ajJ,a'jJ'). (47) 

OLjJ Oi'jJ' 

Here, /i ex (<aj J, cdj J 7 ) indicates the term which only in¬ 
cludes the contribution from the set of the intermediate 
states (oij J, ot’jJ'). The contribution from each such pro¬ 
cess can be measured by the width w of the eigenvalues 
of The widths w for the full exchange 

Hamiltonian H ex and those within l/U approximation 
are shown in Figure [3ja) and (b), respectively. In com¬ 
parison with the contributions to the full Hamiltonian, 
those from the high energy states (A E « 5-10 eV) are 

exaggerated in Hll U . 

In the case of bridge, the main exchange coupling 
arises between the f xyz orbital of Dy and the unpaired 
electron occupying the lowest unoccupied molecular or- 
bital (LUMO) of N|“ (Fig. gd)). The LUMO level in 
Nj - has significantly higher energy compared to the or- 
bital energy of 4/ electrons in Dy 3+ . On this reason and 
also due to a larger space distribution of the LUMO com- 
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TABLE II. K™, and K™,. 


rd) 

'mo 

4.4458 x 10“ 3 

r( 2 ) 

'mi 

2.2245 x 10“ 3 

rd) 

^02 

-3.7903 x 10“ 3 

r ;( 2 ) 

^13 

-3.2216 x 10“ 4 

r;d) 

/Vo4 

1.7181 x 10“ 4 

r ;( 2 ) 

''■'15 

1.5888 x 10“ 5 

rd) 

/v '06 

-1.7242 x 10“ 6 

r ;( 2 ) 

''■'17 

-1.8049 x 10“ 7 

irU) 

'mi 

1.0908 x 10“ 2 

r( 2 ) 

/V 22 

4.8355 x 10“ 5 

rd) 

'm3 

-1.1514 x 10“ 3 

r ;( 2 ) 

/V 24 

-4.6269 x 10“ 6 

rd) 

''■'15 

3.5361 x 10“ 5 

r( 2 ) 

'^26 

7.4785 x 10“ 8 

rd) 

/ m7 

-2.9084 x 10 -7 


4.3815 x 10 -5 

rd) 

'm2 

2.9710 x 10“ 3 

4? 

-2.1423 x 10“ 6 

rd) 

'^24 

-1.3026 x 10" 4 

r (2) 

^37 

2.3716 x 10“ 8 

rd) 

'^26 

1.2614 x 10“ 6 

r ;( 2 ) 

/V44 

4.4502 x 10“ 7 

4V 

1.1745 x 10“ 4 

4V 

-7.2365 x 10“ 9 

4V 

-3.3571 x 10“ 6 

4V 

1.0531 x 10“ 7 

rd) 

''■'37 

2.6060 x 10“ 8 

r ;( 2 ) 

^57 

-1.1732 x 10“ 9 

irU) 

A-44 

5.4544 x 10 -6 

r( 2 ) 

'^66 

1.1839 x 10 _1 ° 

rd) 

''■'46 

-4.9672 x 10 -8 

iA?) 

/V77 

1.3155 x 10 -11 

rd) 

^55 

7.4109 x 10“ 8 



rd) 

''■'57 

-3.9138 x 10 -10 



rd) 

''■'66 

4.1115 x 10“ 10 




pared to the 4/ orbitals, the minimal electron promotion 
energy from to Dy 3+ (E/ 21 ) is expected to be much 
smaller than in the opposite direction (E/ 12 ). Hence, we 
neglect the latter process. Given that the Dy orbitals 
involved in the electron transfer have \m\ = 2, according 
to Eq. (38) (/max = 5, the same for the maximal power 
of J\± in the exchange Hamiltonian. 


Figure [2jf ) shows the exchange levels obtained for full 
i/ ex , its first-rank part, and for the 1/U approximation. 
In the present case, the first-rank part of i/ e x coincides 
with Eq. while the corresponding spectrum strongly 
differs from the full i/ ex indicating the importance of 
higher order terms. As in the previous example, the 1/U 
approximation modifies the relative contributions to the 
exchange interaction from intermediate states (Fig. [ 3 ] 
(c),(d)) and induces, in particular, the interchange of the 
three-fold degenerate ground and the nondegenerate first 
excited states (marked with arrow in Fig. |2jf)) [52] [53] . 
Because of the difference in the nature of the exchange 
states, the magnetic properties predicted by the exchange 
states of the full Hamiltonian and the l/E/ approximation 
differ from each other. 


There is another reason that the 1/U approximation 
is not recommended: the Hund’s rule coupling is com¬ 
pletely neglected within this approximation, leading to 
the removal of the Goodenough’s ferromagnetic exchange 
contribution [40] though the latter plays important role 
in many systems. 


VI. CONCLUSION 

The main results of the present work can be summa¬ 
rized as follows: 

1. We derived the Hamiltonian of exchange interac¬ 
tion between J-multiplets (J-J) and between J- 
multiplet and isotropic spin (J-S) on the basis of 
a complete electronic Hamiltonian, including the 
intrasite relativistic effects. The exchange parame¬ 
ters are expressed via microscopic quantities which 
can be extracted from first principle calculations. 
Despite their microscopic character, the obtained 
expressions ( |33| ) and ( |35| are general (i) for ar¬ 
bitrary choice of quantization axes on two mag¬ 
netic sites (which are not expected to coincide) 
and (ii) for various magnetic ions, which can be 
lanthanides, actinides, transition metal ions under 
special conditions or any their combinations. The 
only requirement is that the low-lying states on 
the sites are well approximated by crystal-field split 
eigenstates of a total angular momentum. 

2. The structure of the J-J and J-S exchange Hamil¬ 
tonian is clarified on the basis of derived exchange 
Hamiltonian. More specific, the maximal rank and 
the projections of the irreducible tensors appearing 
in the exchange Hamiltonian are elucidated. 

3. The obtained form of the (kinetic) exchange Hamil¬ 
tonian was analyzed for different geometries of the 
bridge. The relation between the geometry and the 
structure of the Hamiltonian was established. 

4. On the basis of considered examples, we found that 
the exchange spectrum in systems with J-J and J- 
S interaction cannot be adequately described nei¬ 
ther by exchange Hamiltonian of isotropic form 0 
nor within the 1/U approximation. 

5. The contributions to the kinetic exchange Hamil¬ 
tonian from the intermediate J-multiplets are ana¬ 
lyzed. It is found that the 1/U approximation ex¬ 
aggerates the terms from the excited states. More¬ 
over, within the 1/U approximation, the term split¬ 
ting which is larger than the average U is neglected, 
leading to the wrong order of exchange levels. 

In combination with ab initio and DFT extraction of 
microscopic electronic parameters, the microscopic ex¬ 
change Hamiltonians derived in this work can become a 
powerful tool for the investigation of strongly anisotropic 
materials containing metal ions with unquenched orbital 
momentum. 
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Appendix A: Theoretical tools 


The transformation of the exchange Hamiltonian into 
the tensor form is done using the theory of angular mo¬ 
mentum HHnmEU For the convenience of the readers, 
the tools necessary in the derivation are collected here. 


1. Coupling of angular momenta 


For the phase of the spherical harmonics Y//#, 0), we 
use the convention in Refs. [38] and M With this phase 
convention, the complex conjugation of Y™ is related to 
as 


[Y™{6,ct>)]* ={-l) m Y~ m {e,ct>). (Al) 


Consider two systems whose states are the eigenstates 
of the angular momentum, | jiirii) (i = 1,2). The cou¬ 
pled state characterized by the total angular momentum 
can be constructed using the Clebsch-Gordan coefficients, 

fijm 


I jm) = Yi \hmi,i 2 m2)C]™ mihm2 - ( A2 ) 

mi,m 2 


The Clebsch-Gordan coefficients have following symme¬ 
try properties (Eqs. 8.4.3. (10), (11) in Ref. [38]) : 


fi3sm 3 _ / i \ii+i2-j3/^i3m 3 

jim 1 j 2 m 2 V ) j2m 2 . 


j 2 m 2 j 1 m 1 


n. 


_ ( i \,7i —mi j3 /~ij 2 'm 2 

\ J ] F ^jimijs—m 3 

Ai j2 

_ /_i V7i—mi n,3 y^jj 2 m 2 

— \J tt ^j 3 m 3 ji — mi 


; n 

n 

n' 


^_!ji2+m 2 ll 33 (jji— mi 


3i 


j3 —m 3 j 2 m 2 


_ f i V 72 +m 2 ±J -J3 fyjirrn 

\ J tt . j 2 -m 2 j 3 m 3 

L1 Ji 

= (-l) jl+j2 ~ h C J >I ™ 3 lj2 _ m2 , (A3) 


Here, the subscript j indicates the rank, the superscript 
m is the component, 0 and </> are the spherical angular 
coordinates. 


where = a/2 j + 1. 

Using the Clebsch-Gordan coefficients, the 6j and 9 j 
symbols are defined as [38] 


A . .,A ,( — l)3i+32+h+3T\. 
u jj u mm V ■ L / ■ L± J 


Jl J2 Jl2 

I J3 j J23 


/^iji 2 'mi 2 s-ij'm' ^T/23^23 

^ji2rn 12 j 3 m 3 ^j 1 m 1 j 2 m 2 ^j 1 m 1 j 23 m 23 ^j 2 m 2 j 3 m 3 i V ) 


and 


5 jjf n j 12 j 3 Ajl 3 j 2A 


Jl J2 J12 
j 3 44 434 

413 424 4 


= E 4 

miiriij 


r^3 I2mi2 /^J34m 3 4 m i 3 mi 3 24^24 

Ji2mi2j 3 4m 3 4 Ly iimi42m2 Ly j 3 m 3 4 4 m4 Ly ji 3 mi 3 j24m24 Ly iimij 3 m 3 Ly j2m2j4m4’ 


(A5) 


respectively. Here, 


stands for the summation over all mi and rriij (i,j = 1,2, 3,4). The 6j symbol (A4) is 


symmetric with respect to the permutation of columns and the interchange of the upper and lower components of two 
columns (Eq. 9.4.2. (2) in Ref. [38]) . 

From Eqs. (A4), ( |A5| ), we immediately obtain some formulae involving 6j or 9 j symbol. Multiplying both sides of 
Eq. (A4) by Cj^ 1 j 2a m 2 3 an( ^ summ i n g over j'm', we obtain (Eq. 8.7.3. (12) in Ref. | 


(Oil +32+33+3 TT . . n3 r 

V L ) A Ul 2 j 23 Ly , 7 i 


Jimij23m 23 


I Jl J 2 412 
I 43 4 423 


E fijm fiji2m 12 ^,j 23 

^ ji2rn 12 j 3 rn 3 jirn 1 j 2 rri 2 j 2 rr 


m 23 

J2m 2 j 3 m 3 ’ 


(A6) 


m 2 m 3 mi2 


Similarly, multiplying both sides of Eq. (A4) by C'j 1 7m 12 j 3 m 3 C 'i 1 m 1 j 23 m 23 and summing over we obtain (Eq. 

8.7.3. (12) in Ref. 


_ 1 lii +h +33 +3 rr. . ^ijm 

V ± ) J - L ji2j23 ^ji 2 m 12 j 3 m 3 w i 1 mij2 3 m 23 


jm 


I 41 42 412 
I 43 j 423 


_ \ ^ /^ljl2'mi2 S^lj 23^23 

— Z-^^3i r mij 2 m 2 ^j 2 m 2 j 3 m 3 - 


(A7) 
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Multiplying both sides of Eq. (A5) by Cj 1 ^ ni3:/ - 24m24 and summing over j f m f , we obtain similar formula involving 
five Clebsch-Gordan coefficients: 


tt r^3 rn 

li il2j34jl3j24 ( - y j 1 3mi3j24m24 \ 3S .?34 



/^J34^34 syj 13^13 /^:?24?'™24 f A H) 

Ji2m 1 2j34m34 Ly jim 1 j2m2 U j3m3j4m4 Ly jim 1 j 3 m3 Ly j2m2j4m4' 


Multiplying both sides of Eq. (A5) by and summing over jmj'm ', we obtain a formula 

involving five Clebsch-Gordan coefficients (Eq. 8.7.4. (26) in Ref. 


3 1 J2 J 12 

tt V- ) ■ ■ 

ll ji2j S 4ji3j24 Z^^ji2m 12 j34m 3i ^j 13 m 13 h4m 2 4 S 33 34 J34 

jm [jl3 314 3 


E /^ljl2'mi2 S~lj34'rri34 /^jl3^13 ^1324^24 / A 

^jirn 1 j2rn2^j3rn 3 j4rn4^jirn 1 j 3 rn 3 ^j 2 rn2j4rn4’ V 


2. Irreducible tensor operator 


aIld ^292 fulfillS 


The irreducible tensor operator is defined as the op¬ 
erator T kq which transforms as spherical harmonics Y£ 
(Al) under 50(3) rotations: 


RT kq tf = ]T Tk q 'D q , q (R), 

q' = — k 


(A10) 


where D q/q (R) = (Y£ \R\Y£) is the Wigner O-function 
[38] . R £ 50(3) , and R is the rotational operator for R. 
Since Eq. (A10) holds for any infinitesimal rotations, T kq 
satisfies 

' = y/k{kTY)C k k ^;f kq+ ^ (All) 


Jun T} 


kq 


The matrix element of T kq with respect to the eigenstates 
of the angular momentum, {| JM) }, is proportional to 
the Clebsch-Gordan coefficient ( |A2| ). The Wigner-Eckart 
theorem reads (Eq. 13.1.1. (2) in Ref. | 


(JM'\T kq \JM) = 


{-l) 2k {J\\T k \\J) JM . 

A 


n. 


JMkq- 


(A12) 


One of the irreducible tensor operators is the spherical 
tensor operator Y£( J) which is constructed replacing the 
coordinates in the spherical harmonics Y£(r/r) by the 
total angular momentum operator (r/r J) and aver¬ 
aging it over all possible permutations of J operators [8]. 
For example, is replaced by (J M Ju + Ju J/x)/2. 

When the system consists of two subsystems, double 
tensor issibu is used. In this work, the subsystems are 
the orbital and the spin parts of the system. The or¬ 
bital and spin subsystems transform as irreducible ten¬ 
sor within 50(3) and 51/(2) operations, respectively. 
Thus, for the rotation of the system R = R\R 2 E 
50(3) <g> 517(2), where R x E 50(3) and R 2 E 51/(2), 
the double tensor T^ q ^ of ranks k\ and k 2 transform as 


RT^3R ] = 
&2<?2 


k\ k 2 

E E 

q[ = -k i q ' 2 = -k 2 


f, fel q }D k } (RAD 1 *? (R 2 ), 
& 2<?2 V ' <? 2 <? 2 V 


(A13) 


f rfikiqi 
^dfc 2 g 2 




kiqi 

k 2 Q2 


= y/hfa + l)C k k ^+;f^, (A14) 
= VWh + 1 )C%£tff k k £ +li . (A15) 


One of the double tensor operator is the electron cre¬ 
ation operator in atomic spin-orbital (ra, <j), cj ma . This 
is clear since cj ma creates the one-electron state which 
transforms as the product of the spherical harmonics, 
\lm, \<j). On the other han d, the anni hilati on opera¬ 
tor Ci ma does not fulfill Eqs. 

Ci-m-rj defined below does 


A14) and (A15), whereas 


= (-i) 


Zi+m+|+cr- 

—m—a 


(A16) 


Thus, Ci-m-o- instead of the annihilation operator is a 
double tensor. 


3. Method of equivalent operator 

Consider an irreducible tensor T kq of rank k and its 
argument q acting on the spin degrees of freedom. Re¬ 
placing J in Eq. (A12) with 5, we obtain an expression 

On the other hand, 


of the matrix element of Tkq- On the other hand, the 
matrix element of the spherical tensor operator Y£( S) is 


written as 


(SM'\Yg(S)\SM) = 


(~ir(S\\Y k (S)\\S) SM , 

w 


n.< 


SMkqi 


(A17) 


where S is an abstract spin operator. Comparing Eqs. 
dA12l > and ( |Al7| >, one finds the relation between tensor 
operators T kq and Y)?(S): 


(SM'\T kq \SM) = 


(SlfollS) 

<5||y fe (S)||S) 


(SM'\Y q (S)\SM). 


(A18) 
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This equation holds for any matrix element, and hence, 
in the space of {|5M)}, 


' kq 


(S|+||S> 




(A19) 


(5||r fc (S)||5) 

The reduced matrix element in the denominator is sim¬ 


plified using Eq. ( |A17| ) with M = M' = 5 and q = 0, 

(A20) 


<s||r*(s)||s> = ( S^ Yk(s). 

°SSkO 


Consequently, T^ q is expressed as 

(-l ) 2fc (5||f fc ||g) ss F fc g (S) 
n s SSk0 Y°(S)' 


Tkq = 


Eq. (A21) holds for 


0 < k < 25. 


(A21) 


(A22) 


Y£(S) in the denominator of Eq. (A21) is the scalar 


obtained by substituting S 2 = 5(5 + 1) and S z = 5 in 
the spherical harmonic tensor Yj?( S). 

Now, we consider the case of double tensors , of 
rank (k,k f ). Assuming Eq. (J2|, it is transformed into the 
tensor form within the space of the ground J-multiplet, 
{| JM)}. The matrix element of T*? q , is 


iJM' 


<JM'|T+|JM> = + + C J LML sm'^ LMlSMs 


riJM 

i's Cl 

x {LM' L SM' s \T k ?,\LM l SMs). 


(A23) 


The Wigner-Eckart theorem (A12) is applied to the or¬ 


bital and the spin parts of the double tensor separately, 


(JM'\T k ?AJM) = 


(~1)*W (LS\\ffr\\LS) 


L k'q‘ 


n ls 


E \ " ^ JM' 

/ , °LM^M+LM l 5M s 


-yJM 


is replaced by the matrix element of the irreducible tensor 
operator Y™( J), and , within {| JM)} is expressed as 


= (-l) 2k+2k '(LS\\f k \\LS)J2^J n 

rim 

Yn{ J) 


v rrJJ rrnm 
X ° JJnO^kqk'q' 


L S J 
L S J 
k k' n 


W) 


(A26) 


&, fc', and n in Eq. (A26) obey 


0 < k < 2L, 0 < k' < 25, 

max[|fc — k' |,0] <n< min[fc + fc', 2 J\. (A27) 


In the derivation of Eq. ( |A26 ) we used Eq. Q, while it 
is not mandatory. However, without using Eq. (J2|, the 
result will have a more complicated form. 

In the previous works, the so-called Stevens operator 
o q k ( J) = AkqY£(J) has been used instead of Y£ where 
Ak q is a coefficient which depends on both k and q [42] . 
However, the original Stevens operator does neither obey 
Eq. (All) nor the Wigner-Eckart theorem (A12). In or¬ 


der to use the Wigner-Eckart theorem, one can introduce 
such coefficient A that only depends on k and is indepen¬ 
dent from q [55 , 56 . Furthermore, we write the Stevens 
operator in the form of 0^(J)/0^(J) which is equal to 
Y£(J)/Yj?(J). Since the constant A\~ is canceled in this 
form, it is possible to apply the Wigner-Eckart theorem 
flA12l ): 


< JM,| f+ M> =c§ 


r>JM' 
U JMkq 


(A28) 


JJkO 


In this article, we use the latter (see for example Eq. 
flA26l ». 


Appendix B: Derivation of the exchange 
Hamiltonian in J-representation 


s-yLM'r ^SM'v 
X ^LM L kq^SM s k'q•' 


(A24) 


Using Eq. (A9), the sum of the products of the Clebsch- 


Gordan coefficients reduces to the sum involving 9 j sym¬ 
bol: 

(JM'|T+|JM> = (-l) 2k+2k '(LS\\f k \\LS)^Jn 


v /n iJM' /n mm 

x ° JMnm^kqk'q' 


L S J 
L S J 
k k' n 


>, (A25) 


where n is the rank, m is its argument. The r est proce¬ 
dure is the same as the derivation of Eq. (A21). Cj^ nrn 


1. Direct exchange Hamiltonian 

Detailed calculation of J BE is s hown here. Applying 
the method of equivalent operator (A26), {cj ® in 


Eq. (31) becomes 


{4 ® 5 ;}““ =J2 U k.J,C J r iJ r\ n C kq 


Li Si Ji 

kJi^f i f i k0^aab(3 \ Li Ji 

k( i y a b k 

X (-l) 2a + 26 (L i S i ||{c+c i }j|L i 5 i ) 

. 02(Ji) 


OliJi)' 


(Bl) 
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Introducing V l abk defined by 


V\ 


abk 


n kj i cj\ J J \ ka < 


Li Si Ji 
Li Si Ji 
abk 


X (—l) 2a+2b (I/j5i|| {ct ®Cj}* II LA), (B2) 


Eq. ( |B1[ ) reduces to Eq 

Replace {c\ (g Ci} b a in Eq. (|3l|) with Eq. (32), 




— ( i\Zi+n+i+<7 / \ \ A pikq sp>i 

~ I -1 ! 2^ / ,^aabS U a 

aab(3 kq 


abk 


w Axaa /nfO, 

x u L-mL-n u i 


bf3 


o q K Ji) 




(B3) 


Substituting Eq. (B3) in the direct exchange Hamilto- 

§• 


man 


E 

/eg /c'g' aab/3 a'a'b' /3' 

x ]T y mm ^„(-i) il+n (-i) i2+n ' 

mnm'n' 


era' 

s~ikq piaa pik'q pi a' a' 

A ° aabf3^hmli-n^a'a'b'f3' ° hm'l 2 —n' 


X ^abk^a'b'k’ 


oijWM 

oKJi)ol(j 2 )- 


(B4) 


Since (-1)H CT = (-1)-|" CT , = (-1) _/3 , and 

Eaa' = 8 w 8p,-p>, using Eq. Jill, 


The maximum of b is 1 because Si > 1/2 for the magnetic 
ions. For given (a, 6), k is at most a + b. Simultaneously 
k is less than or equal to 2 Ji (Bl). Therefore, the range 
of k is 

0 < k < min[2/^ + 1, 2 Li + 1, 2«/*]. (B7) 

2. Kinetic exchange Hamiltonian 

In the kinetic exchange Hamiltonian Pke, the op¬ 
erators appear as the form of and 

CimaP^+j 1 c\ m ' a ' - One should note that the projection 
operator P^ j j * s totally symmetric within 50(3) group, 
whereas reducible within the 50(3) (g) 51/(2) group. 
Thus, PP^jj is reduced within 50(3) (g) 51/(2) group in 
order to simultaneously treat it with the other double 
tensors. With the use of Eq. the projection operator 
P&jj = Em \iNajJM)(iNajJM\ is 


oN _ n JM r iJM 

ioijJ — 2_^ 2_^ 2_^ ° LMlSMs^ LM' l SM' s 

M M L ,M s M' l ,M' s 


x |^a LS LM L 5M s )(i^LsEM/5M^|.(B8) 

Introducing the irreducible double tensor PP jaQ , a / Q; / de¬ 
fined by 

piV _ \ A \ A f _-j \LK-m' L -\-S-\-m' s s^iaot 

r ia.jaaa'a' / j / j V L ^LmLL—m' L 


X C a s^ sS _ m , s \Lm L Sms){Lm' L Sm' s \, (B9) 
the projection operator P^jj is written as 

\-m' T +<S'+m / c / 

y LuilSitls 


P Lj = E E E (-1 K+m' L +S+m' sC JM 


dN 


y. piJM \ ^ rraa r^a'a' pl\ 

A ^ Lm' L Sm' s / j ^LrriLL—m' L ^SmsS—m'g^iotjaaa'a' 


HdE = ~ E Y. V k a JAV 2 a'bk' 

kqk'q' aa'b 


2 OK3l)Ol^2) 


o°KJ i)0°vW 


(B5) 


The coefficient of the operators is Eq. ( [33] ) . 

The ranks for the orbital (a) and the spin (6) parts of 
{ct are bounded by 2/^ and 2x1/2, respectively. 

Moreover, from the 9j symbol in Eq. ( |B1[ ), a < 2Li 
and b < 2Si where and Si are the L5-term for the 
ground J-multiplet states, respectively. Thus, the ranges 
of ranks a, b are given as 

0 < a < 2min[^,Li], 0 < b < 2 min[l/2, 5$], (B6) 


(BIO) 

Using the s ymm etry properties of the Clebsch-Gordan 
coefficients (|A3|) and the 6j symbol (A4), 


pN _ V^/_1 \ J-L-S+ayr I 5 5 a I AW 

± iajJ / g *-) p J \ ± ^ ajaaa ~ 

aa l J 

(BH) 

The range of the rank a in Eq. jBlTt is 

0 < a < 2 min[5, L], (B12) 


and —a < a < a . Substituting Eq. jBTTt into 


C'imcfPiajJ ^ina' ( 1 ) 


li +^+ \ +cr / U Ip,. 

^imer 1 iajJ Oi—n—cr 5 


U pNi-1*' _ S^(_i\J-L-S+a tt 


JJ U L J I (- 1 ) , ‘ +n+4+ff, cL.^ <_1 


C P 4 , 

imcr ia,j aaa—a l ~ n ~ <J 


(B13) 
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The operator in Eq. (B13), clmaP^jaaa-a^-n-a', is reduced as follows: 

pNi—i - _ sr n b P nd8 -t f pNi—l ~- \ 

7 r ia jaaa—n— & 7 > ^ aali—n^ a—cx\ —cr' u ima | r iotjaa ^ J 


b(3 

dS 


b/3d8 

= E (B14) 

bPc'ydSee ^ ^ ee 


Using Eqs. (A3) and (A7), 




ima Jr 'iajaaa—a^~ n ~ cr ' 


E C,-„^E(- 1 ) l+2a+ “ +d+/n - 

bfdc^dee fcf) 


\ \ f 
a e d 


C«„c«, {4 0 {p»‘-l ® <%}*}" 


(B15) 


Here, the ranges of the ranks 6, c, d, e, / are 

| a — Z$ | E Z) E ^ T Z^, | n — 1/21 < d < a + 1/2, | b — Z$ | E c E Z? H - Z^ , 

|d — 1/2| E e < d+ 1/2, max[|a — e|, 0] E / E min[a + e, 2 x 1/2]. 

and their arguments satisfy —b < /3 < b , —d < S < d, — c E 7 E c, — e E e E e, — / E 4> E /, respectively. Note that / 
is at the largest 1. Substituting Eq. (B15) into Eq. (B13), 


(B16) 


= (-i)' i+n+H<J, (-i) J_i “ s n JJ E (-i) Hd+/ n d e 


aab(3c'ydeef cf) 


[s S a 

I L L J\ I a 


! 1 4}{4®• (B17) 


Similarly, the operator for the other site in Hke (28) becomes 

Cjm’aP^j,c] n , a , = (-l)h+™' + !+-(_l y'-L'-S'nj'j' £ (-I)"* 


~ h+d n de 


aab(3c'ydeef cf> 


\S’ S' alu \ f 
\L' L' J'(\a e d 


• cz ljn ^_ m , bp cii a c^_ a \cj ® ® an:};;. (bis) 


The operators in Eq. (|B17 ) is written in terms of the total angular momentum usi ng th e method of equivalent 
operator. Applying Eq. (A26) to the irreducible tensor, {c\ 0 {-P^jaa ® Eq. (B17), 


= (-l) U+n+1 *+°\-l) J - L - S Y. (-l) Hd+/ n JJde 


aab/3cjdeef cf) 

\ S S <^ 1)2 2 f l c cl r f(f> V tt t 7 < 

X [L L JJ [a e d J J^o°c 7 ee 

x (—l) 2 c+ 2 e (I/ i S' i || jet 0 0 Aj || J L i 5 i )^|U. 

From Eq. (A27), c, e, and & satisfy additional conditions: 

0 E c E 21/^, OEeE 25*, max[|c — e|, 0] E k < min[c + e, 2 J*]. 

The Clebsch-Gordan coefficients are replaced by the sum involving 9j symbol (A9): 

Pil’jjCinv' = (-l) i, + 2+<T EEE( _1 ) 2/C, hE C, /|7 i m C 'iti-<r'“^« 


T/ Si Ji 
U Si Ji 
c e k 


xt j— - 
imer iocjJ 


oi7i> 


f<f> kq x£ 


ajJfxkQO^j.y 


(B19) 


(B20) 


(B21) 
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Similarly, 


£. , p N o+ 1 z t 

^jm'cr-Lju 


= (B22) 

frk Ur, rr-C 


f<t> kq x£ 


Here, T l aj j Sxk in Eq. (|B2l|) is defined by 


ocjJfxk 


= E(-d 


J—L — *5-1-1 -\-b-\-c-\-d—k 


n J J k kJ t b c cl ef Cj’ j’ Q 


Ji J-S' S’ alJi J / 


abcde 


X (- 


l) 2c+2e (L i S i || |cj ® {^"^0*}*} HLiSi), 



(B23) 


and G l aj j fxk is obtained by replacing the reduced matrix element in Eq. (B23) by (LiSi\\{ci<g>{P^ i ^<g>cl} b d }l\\LiSi). 
The ranges of x and its component £ are 


1/ - h\ < x < f + h, 


(B24) 


and —x < £ < x, resp ective ly. 

Substituting Eqs. (B21) and (B22) into i^KE (28), and first summing over a and cr', the operator part of the 
numerator becomes 


E (cj^p^cip = -EEEEE(- 1 t i ^" /+ 9 '^^i J w 

<j<j' fcf) kq x£ k'q' x'Q 

/~ixi f-ixi fix '£' i' O k ®i)Ok' (^i) fTDoc\ 

^linkq^f^lim^lj—n'k'q '^f ° 


The kinetic exchange Hamiltonian 

^ke = eEEEEE 

i^j kq k'q' fxx' oljJ a' T J' 


IS 


v V V fi fi* n x t r< x Z n x 'i' n x 'i‘ 

2-^mn 2-^m'n' L J L mm' L n'linkq^ f(f)lim^lj—n'k' q'^ f — <f>lj—m' 


x P g j gu^ouh) 

•fcjjjxky^xfx'w o%{Ji)o%, (jj )' 


The numerator is replaced by {t x tYkqk'q' j 36 - aIld We 
obtain Eq. (35). 

The range of k in Eqs. (|B21 ) and ( B22) is k < min[c + 
e,2from Eqs. (B12), (B16), and (B23). Since 
0 < / < 1, c + e < 2a + 21 i + 1, and x <h + 1, the range 
of k becomes 


0 < k < min[2^ + 1, 2 Jj\. 


(B27) 


The range of q is restricted by the transfer parame¬ 
ter as well as the maximal k (fc max ). Considering the 
conservation law for the arguments of Clebsch-Gordan 
coefficients in Eq. (36), q = <j) + m + n, and \<j>\ < 1, 


\q\ E min[^max5 2?7T/ max T l], 


(B28) 


where ra max (> 0) is the maximum projection of the mag¬ 
netic orbital that contributes to the electron transfer. 


Uti + AE^+AE^, 


ja'J' 


(B26) 


Appendix C: Property of J kqk k 


By using the Hermiticity of H (29) and 


t 


we obtain 


ok j) * (-lporm 


(- 1 ) q q (Jkqk'q'T = Jk-qk'-q 


(Cl) 


(C2) 

On the other hand, using the time-reversal symmetry of 
H and 

eoi^e - 1 = [o«(-J)]* = (-i)*-«o*«(j), (C3) 
where 0 is time-reversal operator [8], we obtain 

(- 1 ) fc+fc '-«-«' (J kqk 'q')* = Jk-qk'-q'- (C4) 

Comparing Eq. ( |C2| ) and Eq. ( |C4[ ), 

(-l) fc+fe ' = 1. (C5) 

Therefore, both of k and k' are even or odd. 
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Appendix D: Exchange Hamiltonians for J-multiplet 
interacting with isotropic spin 


The exchange Hamiltonian between J-multiplet and 
isotropic spin is obtained replacing orbital angular mo¬ 
mentum of the spin site (i = 2) with zero. When the spin 
state consists of some nondegenerate molecular orbitals, 
the orbital indices r are introduced. 


1. Direct exchange Hamiltonian 


Since l 2 = 0 and L 2 = 0, J 2 = 52, o! = 0, b = k\ 
and the 9 j symbol in T> 2 , hk , (B2) reduces to 1/n s 2 s 2 k'• 
Therefore, V 2 a , bk , becomes 


V 2 


rr'k 1 


= c 


S 2 S 2 
S 2 S 2 k'0 " 


(-l) 2k '(S 2 \\{ci®c 2r ,} ki lp 2 ) 


n S2 


(Dl) 


Here, a 1 and b are omitted for simplicity, and the molec¬ 
ular orbital index r is introduced. On the other hand, 


'yarr / 

* / kqk'q' 


\h+n-q'jf riaa fikq 

/ > / J V rnrr'n'^l 1 ml 1 -n'^aak'-q'i 

mn a 

(D2) 


where, m! = n' = 0, a' = 0, and q' = —/3 are used. 
The direct exchange parameter between J-multiplet and 
isotropic spin is given by 


/ 7 -DE 
kqk' q' 


(D3) 


From Eq. (B7), the rank for the spin site k! is 0 or 1. 


2. Kinetic exchange Hamiltonian 

Since l 2 = 0 and L 2 = 0, J 2 = 5 2 , a' = b' = d = 0, 
x' = /, e! = fc' and d! = 1/2 from the 6j and the 9j 
symbols in T ( |B23[ ). The values of the 6j and the 9 j 
symbols are (-1) 25 /Tl S ' : (-l) 1+/c 7 n /| j and 


l/S 2 S 2 Jfe', respectively. Thus, ^jjfxk becomes 

T 2 — ( — 1 \ 2S ns S1S2S2 ( _1 

^Oisrr'k’ ~ \ L ) n ^S 2 S 2 k'0\ L ) 

02 

x (S 2 || {4. ® } fe; ||S 2 ). (D4) 

Here, the rank for the orbital part in the double ten¬ 
sor projection operator is removed, hreducible 

tensor form which acts on spin state: 


Q 2 ,, 

27 asrr'k 


P 2 N asaa = E (-l) S+m 4g“ sS _ m ,|5m s )(5m^|. 

ms'm' s 

,(D5) 

,/ fe / is obtained by replacing (5 2 ||{cJ> r 0 
With (5 2 ||{c 2r ®P 2 ^+ 1 4 r ,}^|l^2). On 
the other hand, {f x t}P qk , q , (36| reduces to 

x os'v = e EC4 

mm' ^ 

X ^im'/cg^/c^-g'Zim* (^6) 

Therefore, the kinetic exchange coupling parameter is ob¬ 
tained as 


1/ x t\ xrr ' T 1 G 2 

/7KE v-^ W L Skqk'q , ' / ajJk'xk^a' s rr'k' 

kqk'q' / _j tt A . a 77^2 + 1 

V 0/ Ul 2 + ^lcy r.T + ^ 


xrr' oljJ a' S' 


1 OijJ 

\ xrr' r> 1 


U x A :crr ' T 2 

\ —> v—> kqk'q'^ olj Jk'xk*' a 1 T S' rr'k' 

EE E 


^ V 21 + A e£% + A Egrj, 


(D7) 


From Eq. (B27), the rank for the spin site k' is 0 or 1. 


Appendix E: Reduced matrix elements of the 
creation operators 


In order to calculate the exchange interaction pa ¬ 
ramet ers, the reduced matrix elements in Eqs. (B2), 
(|B23|) must be evaluated. In the direct exchange 


interaction (33), there appear (— l) 2a+2h {LiSi\\{c\ 
CiVjLiSi) and (-l) 2a + 2b (L j S j \\{c j ®cl}$\\L j S j ). How- 
ever, note that the irreducible tensor operators are 
the same type as V 11 used for the calculations of 
the spin-orbit coupling m- On the other hand, for 
the calculations of the kinetic exchange interactions 

t ^ f £) N i-l ** 7i.\b 1 


Q, (—1) 2c+2e (T j S'* 11 {cj ® {P^-^c^yjLA) and 
(-l) 2c+2e {L j S j \\{c j ®{P?J+ 1 a ®c]} b d } c e \\L j S j ) have to be 
evaluated. By straightforward calculations, the former is 
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TABLE III. Reduced matrix elements of the creation operator for the ground 6 if-term of Dy' 


.3 + 


r 

(-1)W+W 

|/ 8 ,r)/n iS r 

(-l) 2(l+1/2} (f 9 , 6 H \cl\f 8 ,T)/U L s 

r 

(- 1 )2( i+ l/2)( / 10, r |ct|/ 9 , 6 J7>/ rIis 

7 f 

Vl 

5 G(2) 

1 / 65 

7 V 11 

5 d 

_2 ,/IT 

Z V 35 

5 D(1) 

1 /10 

3 V 7 

5 G(3) 

3 [3 

7 V 2 

5 f 

fM 

V 35 

s D{ 2) 

2y/b 

21 

5 H( 2) 

-1 

5 g 

fl3 

V 7 

S D( 3) 

3 [lb 

7 V 11 

5 J(1) 

1 /FT 

3 V 11 

5 / 

fl4 

V 5 

5 F(1) 

1 

V6 

5 7(2) 

1 f~26 

3 V 11 



5 F(2) 

[Y 

V 22 

5 K 

[lb 

V 11 



5 G(1) 

/ 65 

V 154 

5 l 

fl7 

V 11 




Hi; Si 

(_l)2(G + l/2 )yNi 





^LiMLCj^SiMsee 


I T-LiSi 


v \ \ \ ^ (_-\\li+n+L+NL r^ c l rhfi ^iaa ^iLN l 

X V L ' ° l i mb^ao l .li-n^LN f L L-N L ^L i M , L l i m^L i M L l i 

mn af3 N L N' T 


xEEE(-D 

crcr' a'6 NsN' s 


b -\-S-\-N s rrec 


C ee /o idS ri 


sySN s ^ySNs 

bvdS^aa'b-cr'^SN^S-Ns^SiM'lm^SiMsbn’ 


(El) 


and the latter is 


11 LjSj 


^Lj M' l ^Sj M' s 
^LjMLC^SjMsee 


(_l)2fe + l/2) {l N j+ 1 5 QL 5|| g t ||^ , £ .£.) 


n_L5 


v \ \ \ ^ (_-\\lj+m+L+NL syc'y s^b/3 >oaa ^<LN l 

X ' ^Ij-mb^aodjn^LN^L-NL^LjM'Jjm^LjMLljn 

mn ctf3 N l N' l 

,SN' s 

SjMg 2 / 


v \ A ^ \ A ( -\\b+°+S+Nsr^dd r>iaa' rr SN s r<SN s 

X 2^2^ V ^ ^h-vdd^aa'^'^SN^S-Ns^SiM^mSjMsW 


crcr' a'S NsN' q 


(E2) 


Here, the components Ml, M' l , Ms , M^, 7 , e are chosen 

so that ^L*Mfc 7 CsImscc / 0 is satisfied. For the cal¬ 
culations of the equations, the reduced matrix elements 
of the creation operators, cj ma , are necessary. The are 
calculated as El 


where ( f N aLS{\f N ~ 1 (a'L'S')fLS) is the coefficient of 
fractional parentage, which is tabulated in Ref. 08l The 
reduced matrix elements of c* necessary for the present 


examples are shown in Table III 


Appendix F: Exchange states of Nj bridged Dy 
dimer 


.3+ 


(_l)2(i + l/2)^iy aZ; g|| g t||/iV— Iq/j/g') 

IIls 

= v r N{f N aLS{\f N -\a'L'S , )fLS), (E3) 


The difference between the eigenstates of the full ex¬ 
change Hamiltonian lT ex , |\fq, T), and those of the Heisen¬ 
berg type Hamiltonian (HI), | J 12 M 12 ), are compared as in 
the case of the linear system. The low-energy states are, 
in the basis of {\J 12 M 12 )} ([46]) , 
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|^i, b lu ) « 0.767|1,0) + 0.493|3,0) + 0.108(|3, -2) + |3,2)) + 0.345|5,0) + 0.118|7,0), 

|^2, a g ) « 0.552|0,0) + 0.661|2,0) + 0.400|4,0) + 0.250|6.0), 

I *3, h g ) » 0.586(— 12, -1) + |2,1)) - 0.260(—14, -1) + |4, 1)) + 0.256(-|4, -3) + |4,3)) + 0.108(-|8, -1) + |8,1)), 
1^4, h u ) « —0.545(—11, -1) + |1,1)) - 0.158(—13, -1) + |3,1)) - 0.161(-|3, -3) + |3,3)) + 0.293(-|5, -1} + |5,1)) 

- 0.218(— 15, -3) + |5,3)) - 0.129(—17, -1) + |7,1)), 

|^ 5 , bi u ) « —0.285|1,0) + 0.680|3,0) - 0.448(|3, -2) + |3,2)), 

l^e, a g ) « 0.446|0,0) -0.331|2,0) +0.347(|2,-2) + |2,2)) - 0.460|4,0) + 0.255(|4, -2) + |4,2)) + 0.181|6,0) 

— 0.155(|6, —2) + |6,2)). 

Here, the irreducible representation F of D 2 h symmetry is used. The exchange states belong to the eigenvalues 
Ex = -0.274979, E 2 = -0.274977, E 3 = -0.262325, E A = -0.262299, E 5 = -0.260815, E 6 = -0.260805, respectively. 
As in the axial system, the low-energy states are not well described by Hueis- Within the 1/U approximation, the 
low-energy states become 

\V{ /U ,a g ) = 0.509|0,0) + 0.696|2,0) + 0.434|4,0) + 0.222|6,0), 

| y\ /U ,b lu ) = 0.740|1,0) +0.559|3,0) +0.330|5,0) + 0.115|7,0), 

\*\ /U , b 3u ) = 0.554(—11, -1) + |1,1)) + 0.205(—13, -1) + |3,1)) + 0.193(-|3, -3) + |3,3)) - 0.234(-|5, -1) + |5,1)) 
+ 0.220(—15, —3) + 15,3)), 

b 2g ) = 0.608(—12, -1) + 12,1)) - 0.188(—14, -1) + |4,1)) + 0.279(-|4, -3) + |4,3)), 

| ^l /U ,a g ) = 0.489|0,0) — 0.263|2,0) — 0.374(—12, —2) + |2,2)) - 0.430|4,0) - 0.265(-|4,-2) + |4,2)) + 0.141|6,0) 

+ 0.127(—16, -2) + |6,2>) - 0.103(—16, -4) + |6,4>), 

\t>l /U ,bi u ) = 0.351|1,0) -0.618|3,0) +0.471(|3,-2) + |3,2)) + 0.109(|5,-4) + |5,4)). 


The exchange states \^ l J U ) are also ordered as the in- the 1/U approximation. |\l/j) and T, 1 3 4 5 6 ^ 7 8 9 10 11 12 ) with the same 

crease of the energy. Some levels are interchanged due to representation quantitatively differ from each other. 
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